Compaction dynamics of a magnetized powder 
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We have investigated experimentally the influence of a magnetic interaction between the grains 
on the compaction dynamics of a granular pile submitted to a series of taps. The granular material 
used to perform this study is a mixture of metallic and glass grains. The packing is immersed in 
an homogeneous external magnetic field. The magnetic field induces an interaction between the 
metallic grains that constitutes the tunable cohesion. The compaction characteristic time and the 
asymptotic packing fraction have been measured as a function of the Bond number which is the 
ratio between the cohesive magnetic force and the grain weight. These measurements have been 
performed for different fractions of metallic beads in the pile. When the pile is only made of metallic 
grains, the characteristic compaction time increases as the square root of the Bond number. While 
the asymptotic packing fraction decreases as the inverse of the Bond number. For mixtures, when 
the fraction of magnetized grains in the pile is increased, the characteristic time increases while 
the asymptotic packing fraction decreases. A simple mesoscopic model based on the formation of 
granular chains along the magnetic field direction is proposed to explain the observed macroscopic 
properties of the packings. 

PACS numbers: 81.05.Rm,81.20.Ev 



I. INTRODUCTION 

Over the last decade, the compaction of granular ma- 
terial has been the subject of numerous studies in the 
physics community [l], 0, El] as the problem concerns a 
large panel of both industrial and fundamental research. 
The majority of these studies are focused on cohcsive- 
lcss granular materials. However, a better fundamental 
knowledge of cohesive powder properties is required for 
their manipulation, in particular within the recent de- 
velopment of nanopowder technology. Indeed, cohesive 
forces are known to strongly affect the flow and the static 
properties of fine powders, since they induce the forma- 
tion of large aggregates [1] . 

The compaction dynamics of cohesiveless granular ma- 
terial submitted to a series of taps has been studied ex- 
perimentally, theoretically and numerically @, @, HBH 
[Tol . ITT1 . [l2l , fl3l ] . Different laws have been proposed for the 
evolution of the volume fraction 77 as a function of the 
number n of taps. Among others, one has proposed the 
inverse logarithmic law [5J 



A?] 



1 + B ln(l + 2) 



(1) 



where the parameters 7700 and A77 are respectively the 
asymptotic volume fraction and the maximum variation 
of the volume fraction. The dimensionless parameter B 
depends on the acceleration during each tap and r is the 
characteristic relaxation time for the grain reorganization 
process. 

The difference between a powder and a granular assem- 
bly is the increase of the importance of cohesion forces 
in respect to the gravitation because the grains are much 
smaller in powder than in granular assembly. It is diffi- 
cult to imagine a model powder for which it is possible to 
tune the size or the shape of the grains independently to 



the cohesion. In order to study the influence of the cohe- 
sive forces on the properties of granular assemblies, sev- 
eral methods can be used. For example, one may change 
the size of the particles [lH . However, controlling the size 
and the shape of particles is a difficult task in particu- 
lar for small cohesive objects. Liquid bridges can also be 
used to increase the interparticle forces 
certain conditions, to lubricate contacts 
to cohesive forces. 
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FIG. 1: (a) Interactions between two ferromagnetic particles 
in a magnetic field B. The /2 vector denotes the magnetic 
dipole induced by the magnetic field B, rtj denotes the vector 
between the dipoles. (b) The magnetization of two ferromag- 
netic particles in a vertical magnetic field induces repulsive 
interaction when r7j is perpendicular to B and attractive in- 
teraction when rtj is parallel to B. 



The idea of our study is the use the magnetic force to 
measure the impact of the cohesion on the compaction 
properties of a granular assembly. The interparticle force 
can be tuned continuously by varying the strength of the 
external magnetic field B. This technique has previously 
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been used to study the influence of the cohesion on the 
packing fraction [l?], [3 and on the avalanche angle [l9[ 
of granular assemblies. The choice of the size of grains 
(beads) is important: not too small to avoid 'natural' 
cohesion, not too large to avoid the use of high magnetic 
field and to remain in powder conditions. 

The influence of the magnetic field can be sum up as 
following. Each ferromagnetic particle becomes a mag- 
netic dipole characterized by a magnetic moment pi. We 
assume that the magnetic moments fl arc parallel to the 
applied magnetic field B. This is an approximation when 
we consider many interacting particles. The potential en- 
ergy Uij between two magnetic dipoles i and j separated 
by a distance (see Figure [T]) is given by 

TT _ mm^ -i - 3cos 2 /3 



where (3 denotes the angle between the vector fij and 
the magnetic field B. As illustrated in Figure [1] the 
potential is attractive when (3 < 55° and repulsive when 
55° < j3 < 90°. These interactions between the particles 
change significantly the internal structure of the packing 
and the force network. Moreover, due to the directional 
nature of the force acting between dipoles, anisotropic 
particle aggregates form along the field direction. 

In this paper, we study the influence of the magnetic 
interaction between the grains on the compaction dynam- 
ics. The size of the grains has been chosen such as the 
cohesion is weak regards their weights and small enough 
to be considered as a powder. First, a granular material 
made of metallic beads has been used. Afterward, mix- 
tures of glass and metallic beads have been considered. 
The compaction characteristic time r and the asymp- 
totic packing fraction rjeo have been measured for differ- 
ent magnetic field strength and for different fraction <fi of 
metallic beads in the mixture. At the end of the paper, a 
simple mesoscopic model is proposed in order to describe 
the properties of a magnetic packing. 

II. EXPERIMENTAL SETUP 

A sketch of our experimental setup is illustrated in Fig- 
ure [2j A glass tube of internal diameter D = 21 mm is 
placed into a vertical magnetic coil where a constant cur- 
rent can be injected. The magnetic field is thus parallel 
to the gravity field. The strength of the magnetic field 
can be fixed between and 150 Gauss. We checked that 
the variation of the magnetic field along the vertical axis 
of the coil does not exceed 10 %. The average diame- 
ter of both metallic and glass grains is rf=100 /im with a 
polydispersity of about 50%. 

In order to obtain a reproducible and spatially homoge- 
neous initial packing fraction, the following initialization 
protocol has been used. A narrower and bottomless glass 
tube (external diameter slightly below 21 mm) is inserted 
into the main vessel. Afterward, the small tube is filled 
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FIG. 2: Sketch of the experimental setup, (a) An initializa- 
tion protocol is used to obtain a reproducible and homogenous 
initial packing fraction 770. A smaller bottomless tube is in- 
serted into the main glass cylinder. This tube is filled with 
the granular material and is removed at constant speed v. 
(b) After the initialization process, the magnetic field B is 
set. An electromagnetic hammer produces periodic taps on 
the bottom of the system. A CCD camera take a picture of 
the top of the pile after each tap. In order to keep the gran- 
ular/air interface flat, a light cylinder is placed on the top of 
the pile. 



with the granular material. The small tube is then re- 
moved upward at a low and constant velocity v = lmm/s. 
The current in the coil is set after the initialization proto- 
col in order to obtain a constant initial packing fraction 
7]q. A light disk is gently placed on the top of the pile in 
order to keep it flat during the compaction process. The 
system used to produce the taps under the tube is an 
electromagnetic hammer as described in || |22j ■ During 
one tap, the system undergoes a short peak of accelera- 
tion (the width of the peak is 0.25 ms and the maximum 
intensity is 15g) and some damped oscillations during 2 
ms. To obtain a compaction curve, two thousand taps 
are produce. Successive taps are separated by one sec- 
ond. Therefore, the packing relaxes after each tap. 

The packing fraction of the pile is measured by image 
analysis from a CCD camera. The average position h of 
the upper interface of the packing is measured after each 
tap and gives the packing fraction r) by using the known 
value for the volumic mass p of the beads and the total 
mass m of grains. One has r\ = -^js^ ■ The part of the 
pile situated above the coil is typically less than 10 % of 
the total height of the pile. 

To obtain a relationship between the magnetic field 
strength and the intcr-grain interactions, a reference 
magnetic field B rc j has been measured. This refer- 
ence magnetic field gives the balance between the mag- 
netic force acting on two contacting grains aligned with 
the magnetic field and the weight of one grain. Then, 
we can compute the Bond number Bo which is the ra- 
tio between the cohesive force F c and the grain weight 
rag. If we assume that the inter-grains force is propor- 
tional to the square of the magnetic field [HI , we obtain 
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FIG. 3: Typical compaction curves corresponding to different 
Bond numbers Bo. From the top curve to the bottom one, 
we have: Bo = 0.0, Bo = 0.4, Bo = 1.0 and Bo = 5.8. The 
fraction of metallic grains is (f> = 1. The experimental data 
are fitted with the logarithmic law EqfTJ. 

Bo = B 2 /B^ e j:. The reference magnetic field for the 
metallic beads is B re f = 50G. 



III. THE PURE CASE: ONLY METALLIC 
GRAINS 

A. Observations 

To perform the compaction experiments, the pile is ini- 
tialized with the method explained previously. The mag- 
netic field is set after the initialization of the pile in order 
to obtain always the same initial packing fraction tjq. The 
figure [3] presents four typical compaction curves obtained 
with different values of the magnetic field strength. The 
corresponding Bond numbers arc Bo = 0, Bo = 0.4, 
Bo = 1.0 and Bo = 5.8. The continuous curves are the 
fits with the logarithmic law Eq. (fTJ). The parameter B is 
fixed to 1 in order to minimize the number of free fitting 
parameters. Indeed, the parameters B and r are known 
to be strongly anticorrelated. 

The evolution of the characteristic time r obtained 
from the fits as a function of the Bond number is pre- 
sented in Figure 2] The strong increase of r observed 
for a low Bond number shows that the compaction dy- 
namics is strongly influenced by the interactions between 
the grains. The cohesion induced by the magnetic field 
slows down the compaction process. We have observed 
previously this effect with fine cohesive powder [22| . The 
evolution r as a function of the Bond number Bo is well 
fitted by a square root law (see Figure[4ja)). Therefore, t 
increases linearly with the magnetic field strength B. In 
the cohesionless case (Bo = and B = 0), some fluctua- 
tions are observed in the compaction curve after 80 taps 
(see figure [3j) . These fluctuations are related to a con- 
vective process inside the packing. When the Bond num- 
ber increases, this convective process takes place later. 
Moreover, for high Bond numbers (Bo > 1), i.e. when 



the cohesion between the grains is higher than the weight 
of one grain, this convection process disappears. 

Figure 0Jb) presents the asymptotic packing fraction 
r/oa as a function of the Bond number Bo. The frac- 
tion rioo decreases when the magnetic interaction between 
the grains strengthens. This behavior has been observed 
in a previous work [l7| for the packing fraction t]o of a 
metallic grains pile created in a vertical magnetic field. 
One should note that the asymptotic packing fraction 
obtained in the cohesiveless case is far lower than the 
random close packing limit t\rcp ~ 0.64. The low value 
of the packing fractions obtained in our experiment is 
due to the roughness of the metallic grains. 



B. Modeling 

In order to estimate the compaction characteristic time 
t, numerous theoretical models [H,[2(| have been devel- 
opped. The great majority of them considers particles 
which are moving in "cages" constituted of neighboring 
particles. As shown recently [2rj| in numerical simula- 
tions, the probability that a particle escapes a cage is re- 
lated to r. The escape probability is a complex function 
of both packing fraction r\ and some energy barrier Ef, 
related to the work for se par ating two grains. Roughly, 
one has a scaling r ~ Ef, [2y|. If one considers that the 
energy barrier is herein linked to the work for separating 
two magnetic dipoles, one obtains the scaling r ~ Bo 
which is not observed in Figure |4j As a consequence, 
the assumption of caging dynamics does not hold in our 
system. 

Theoretical and numerical works have been done to 
predict the packing fraction r\ of an assembly of particles. 
Thermodynamic approaches [23[ and statistical methods 
[24| have been developped. For identical spherical grains, 
an average over all possible Voronoi cells around parti- 
cles allows to estimate the random loose packing frac- 
tion in 2d. A similar approach (2f| allows to estimate 
the packing fraction of bidispcrsc 2d grains. Herein, the 
situation is more complex because the attractive and re- 
pulsive interactions should be taken into account. Note 
that vertical chains are favored by the magnetic field. 

Let us propose a simple model to explain the scaling for 
both compaction parameters r and 7700 ■ Due to the orien- 
tational nature of the interaction, the magnetic field in- 
duces the formation of vertical chains inside the packing. 
The vertical length of a chain \\\ ~ Bo is the key physi- 
cal ingredient for describing the local grains organization 
in the packing (see Figure [5]) . The geometry of a chain 
looks like a disordered aggregate (similarily to Id random 
walk) rather than a linear needle. The "thickness" x± of 
random aggregates is given by x± ~ y/X\\ ~ "^Bo as for 
Id random walk. The packing of such random objects 
involves a loose packing fraction t/q. For increasing the 
packing fraction 77, the structure of the disordered chains 
should change in order to align the grains. Indeed, the 
space occupied by a chain is strongly decreased. Even if 
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FIG. 4: (a) Evolution of the characteristic time r as a function of the Bond number Bo. The fraction of metallic grains is 
0=1. Each point corresponds to an average over three experiments. The continuous line is a fit with a square root law 
t — to + aV Bo. The only fitting parameters is a = 3.8 ± 0.2. The value of the characteristic time when Bo = is to = 10.3. 
(b) Decrease of the asymptotic packing fraction r/oo with the Bond number Bo. The fraction of metallic grains is = 1. Each 
point corresponds to an average over three experiments. Error bars are indicated. The continuous line is a fit by Eq. JSJ. The 
only fitting parameter is 77oo,Bo=oo = 0.388 ± 0.001. 



the height of the chain may increase, several chains can 
arrange in the space occupied by a single chain before 
the compaction process. The characteristic time needed 
to align chains along the magnetic field, being the charac- 
teristic time t for the packing reorganisation, is propor- 
tional to x± > giving the right scaling r ~ V Bo obtained 
in Figure 3J Adding the appropriate limit conditions, one 
writes 



mesoscopic level. Indeed, both size and shape of particle 
aggregates play a major role in such system. 



IV. MIXTURE CASES 



T = T 



(3) 



Observations 



The experimental data are well fitted by the above square 
root law (see Figure S]) with the only one fitting param- 
eter a = 3.8 ±0.1. 

If one consider that the packing is an assembly of cells 
illustrated in Figure [5j the packing fraction r\ of the pile 
is the ratio between the volume of the grains in the cell 

a. One obtain 



Vgrains and the volume of the cell V c 



'/ 



ANttcP 



(4) 



where N is the number of grains in the chain and d is 
the grains diameter. Since N ~ Bo, x± ~ V Bo and 
X\\ ~ Bo, we get the scaling 77 ~ Bo^ 1 . Adding the limit 
conditions, one writes 



7] 00, Bo— 00 



1 + Bo ' 



(5) 



The decrease of rjao with Bo is well fitted by this relation 
(see Figure ). The only fitting parameter is the asymp- 
totic packing fraction obtained for infinite Bond number 
Voc,Bo=oo = 0.388 ±0.001. 

In summary, we have derived a simple model for de- 
scribing the compaction dynamics of a magnetized pack- 
ing. It appears that the packing fraction as well as com- 
paction dynamics are understood when considering the 



The characteristic compaction time r and the asymp- 
totic packing fraction has been measured for mix- 
tures of glass and metallic beads. Eight volume fractions 
<p of metallic beads in the mixture have been considered. 
Moreover, the measurements have been repeated for dif- 
ferent values of the Bond number Bo. 

The figure [5] presents the evolution of the characteristic 
time t as a function of the fraction <p of metallic beads 
in the mixture for three values of the Bond number Bo. 
For each value of </>, the characteristic time To obtained 
when Bo = has been subtracted. As expected, the com- 
paction characteristic time t increases with the fraction <p 
of metallic grains in the pile. Indeed, the cohesion inside 
the packing is expected to strengthen when the fraction 
of magnetized beads increases. This effect is amplified 
when the Bond number Bo increases. 

The evolution of the asymptotic packing fraction rjao 
as a function of the fraction <fi of metallic beads in the 
mixture for different values of the Bond number Bo is 
presented in Figure [7J The asymptotic packing fraction 
?7oo,o obtained when Bo = has been subtracted for each 
value of 4>- The packing fraction is found to decreases 
with the fraction (f> of metallic beads. However, when 
Bo > 1, a local maximum is observed for <j> ~ 0.2. The 
signification of this maximum will be discussed bellow. 
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FIG. 5: (top) An aggregate composed by several metallic 
beads before and after the compaction process. The disor- 
dered nature of the chain formation defines two typical length 
X|| and x± which are related to each other (see text), (bot- 
tom) When the fraction <f> < 1, the chains of magnetized 
grains are divided due to the presence of glass beads. 



FIG. 6: Evolution of the characteristic time r as a function 
of the fraction of metallic beads in the mixture for different 
values of the Bond number Bo. The characteristic time to 
obtained when Bo = has been substracted. Each point 
corresponds to an average over three experiments. Error bars 
are indicated. 
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FIG. 7: Evolution of the asymptotic packing fraction 7700 as 
a function of the fraction of metallic beads in the mixture 
for different values of the Bond number Bo. The asymptotic 
packing fraction 7700,0 obtained when Bo = has been sub- 
stracted. Each point corresponds to an average over three 
experiments. Error bars are indicated. 



B. Modeling 



The previous model can be generalized to mixtures. 
The presence of glass beads in the packing decreases the 
length of the magnetic chains inside the packing (see Fig- 
ure [5]) . The decrease of the chain length as a function of 
the fraction <p can be evaluated simply. If we consider a 
reservoir of metallic and glass balls with a fraction <fr of 
metallic balls, the probability to pick randomly a metallic 
ball is <f>. The probability to obtain n successive metallic 
balls is <j> n . In our system, when (f> — 1, the chains have 



< N >- 



i0' 



(6) 



Therefore, both Equations ([3]) and can be gener- 
alized in the case of magnetic and non-magnetic grains 
mixing. Indeed, in the case of mixtures, the character- 
istic size x\\ is proportional to X^i"" i4> % Bo. Then, we 
obtain 



to 



E 



Bo 



(7) 



average size < N > is 



and 
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FIG. 8: Scaling of the data presented in Figure [6] The plain 
squares, circles and triangles correspond respectively to Bo = 
0.25, Bo = 1 and Bo — 5.8. The data collapse on a same 
trends and are well fitted by the Equation ©. 



T]oo — ^]cc.Bo— oo 



Voo,Bo=0 — Vco,Bo=oa 



-Bo 



(8) 



The division by the arithmetic series ^2f = ^ L ax i has been 
used to recover the Equations and ([5]) when 0=1. 
Both series could be evaluated analytically. 

The Figure [5] shows the scaling (r — T )/\f~Bo as a 
function of <j). The points corresponding to different val- 
ues of the Bond numbers Bo collapse on a same trends. 
This collapse confirms the relation r ~ V Bo. More- 
over, the data are well fitted by the Equation ([7|) with 
a = 4.55±0.7 and N max — 3±1 as free fitting parameters. 

The deviation between the experimental values of 
7700 and the prediction of the model r^oo, mode/ obtained 
from the Eq. ([5]) is presented in Figure [5] The ratio 
Voo/Voo, model is close to unity except for the factions <ft 
of metallic grains situated between 0.1 and 0.3 where a 
local maximum of rjoo is observed. Within this range, the 
metallic grains are mainly in contact with glass grains if 
the mixture is homogenous. However, during the com- 
paction process, some clusters of metallic grains could 
be created in the packing. The existence of this segrega- 
tion process could explain this increases of the asymptotic 
packing fraction. 

In summary, the simple model proposed in the pure 
case has been generalized to the mixture case. We expect 
that the magnetic chains are divided due to the presence 
of glass beads in the packing. 



V. CONCLUSIONS 



The characteristic compaction time r and the asymp- 
totic packing fraction r/aa has been measured for mixtures 
of glass and metallic beads in a magnetic field. When the 
pile is made of metallic grains, the characteristic com- 
paction time r increase as \J Bo and the asymptotic pack- 
ing fraction t]^ is proportional to Bo -1 . We have derived 
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FIG. 9: Deviation between the experimental values of 7]oo 
and the prediction of the model (Eq. ©)■ The plain squares, 
circles and triangles correspond respectively to Bo = 0.25, 
Bo = 1 and Bo = 5.8. 



a simple mesoscopic models that explains these scalings. 
This model has been generalized to mixtures. 

These fundamental results should be useful to interpret 
the measurements made with cohesive powders used in 
the industry. Indeed, the tapping measurement is a very 
common tool to characterize the physical properties of 
powders. 
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